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We analyze the strength of polarization correlations between two light beams that can be achieved 
in the semiclassical regime using statistical mixtures of coherent states and binary on/off detec- 
tors. Under certain symmetry assumptions, the visibility of polarization correlations is shown to 
be bounded by |, which is in a striking contrast with perfect 100% correlations exhibited by pho- 
ton pairs prepared in the singlet state. The semiclassical limit is demonstrated in a measurement 
performed on a pair of laser beams undergoing correlated depolarization. This result illustrates 
the dramatic difference between predictions of quantum mechanics and the semiclassical theory of 
electromagnetic radiation for the polarization degree of freedom. 

PACS numbers: 



I. INTRODUCTION 



Quantum mechanics allows for correlations that can- 
not be understood within the classical paradigm of lo- 
cal realism, which naturally grows out of experiences 
gathered in the macroscopic world Models based on 
this paradigm are not able to reproduce fully predictions 
of the quantum theory unless supplied with additional 
controversial assumptions like non-localityjSuperluminal 
communication, or negative probabilities 0, y, 0, H| • A 
great deal of effort has been dedicated to identify the 
borders defined by local realistic theories, beyond which 
quantum mechanics offers the most dramatic manifesta- 
tions of its unique character. This work currently finds 
applications in novel, quantum-based protocols for secure 
communication and information processing [1, 0, Q ■ 

One of the major sources of excitement in quantum op- 
tics has been tracing effects that can be explained only on 
the grounds of the fully quantum theory of light and mat- 
ter. A standard reference point here is the semiclassical 
theory of photodetection, in which a quantized material 
system interacts with light treated classically. This the- 
ory suffices to explain a wealth of optical observations in 
systems comprising stochastic classical fields, linear op- 
tical elements, and typical photodetectors [9|. A direct 
demonstration of the quantum nature of light requires 
one to resort to more exotic ways to generate and manip- 
ulate light 10]. The semiclassical regime can be viewed 
as a special case of a local realistic theory, and its bounds 
are typically tighter, as discussed in detail for two-photon 
interference effects [2]. However, the semiclassical regime 
can be easily implemented in an optical laboratory using 
readily available resources, unlike other models consid- 
ered in general local realistic theories which may include 
highly speculative assumptions. 

When identifying the limits of the semiclassical regime, 
it is crucial to take into account a realistic description of 
the employed components. For example, most detectors 
with single photon sensitivity do not have photon number 



resolution. As a result, intensity correlations typical for 
non-classical sources of optical radiation can be mimicked 
by suitably engineered non-stationary classical fields, and 
additional tests of quantumness are needed [ll| . A thor- 
ough study of the actual construction details of standard 
photon counting modules can even provide a complete 
control over their behaviour using classical light 
which leads to new eavesdropping strategies for quan- 
tum cryptography [l3l ]. This clearly shows that under- 
standing the precise borders between semiclassical, local 
realistic, and quantum regimes remains full of challenges 
that grow out of previous advances in quantum optics 
and are highly relevant to emerging quantum-enhanced 
technologies. 

The purpose of this paper is to explore the question 
to what degree the semiclassical regime can reproduce 
polarization correlations. Such correlations are a cele- 
brated signature of quantumness [14| and nowadays can 
be routinely observed for photon pairs produced in a 
suitably arranged process of parametric down-conversion 
[15j . Typical for most quantum optical experiments, we 
will assume here that the photodetectors provide only a 
binary response telling whether at least one photon has 
been detected or none at all [HI [l?], [H, Gil • The prob- 
ability of a click is a non-linear function of the incident 
intensity. This can be viewed as a generalization of the 
Malus law Q considered within the semiclassical theory. 
The theoretical discussion is illustrated with a photon 
counting measurement performed on a pair of collectively 
depolarized light beams that provide the strongest form 
of polarization correlations allowed in the semiclassical 
regime. 

The quest to understand the exact relation between the 
classical and the quantum has been vigorously pursued 
by and will remain an important element of the scientific 
legacy of late Krzysztof Wodkiewicz. This contribution 
is dedicated to his memory. 
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II. QUANTUM CORRELATIONS 

The notion of correlations in classical physics is re- 
served to the statistical description of systems for which 
complete knowledge is not available. The state of a sys- 
tem is said to be correlated whenever the joint probability 
distribution describing parameters of the entire system 
does not factorize into marginal probability distributions 
for individual subsystems. If the complete knowledge 
about the system is available, the probability distribution 
is clearly a product of marginal probability distributions, 
since the state of each subsystems is definite and hence 
no correlations are present. 

According to the standard interpretation of quantum 
mechanics, the complete knowledge of a quantum sys- 
tem is encoded in a vector in an appropriate Hilbert 
space. The superposition principle together with the ten- 
sor product structure of composite states implies the ex- 
istence of entangled states — a unique quantum concept 
comprising the complete knowledge about the composite 
system and correlations within. The strength of corre- 
lations present in entangled states may well outperform 
any conceivable classical ones and therefore they provide 
an argument for the incompatibility of local realistic the- 
ories with the quantum world via Bell's inequalities [20I ]. 

One of the simplest and at the same time the most 
characteristic examples of an entangled state is the sin- 
glet state. When realized in the polarization degree of 
freedom of two photons, the singlet state reads: 
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where the kets |<->) and ||) correspond to the horizontal 
and the vertical polarization states respectively, and the 
indices a and b refer to the two photons that are assumed 
to be distinguishable by another degree of freedom, for 
example the direction of propagation. A feature that dis- 
tinguishes the singlet state from other two-photon entan- 
gled states is its invariance under arbitrary polarization 
transformations, provided that the same transformation 
is applied to both the photons. Consequently, the an- 
ticorrelations present in the horizontal-vertical basis in 
Eq. {T]) manifest themselves also in any other basis the 
polarization is being measured. 

More generally, let u a and be two Bloch vectors of 
unit length representing the polarization basis in which 
the polarizations are being measured on the two photons 
respectively. It is easy to verify that the density operator 
corresponding to the singlet state can be written as a 
special case of a two-qubit Werner state [2l| 



g=-\i®i+ri 



i—x,y,z 



(2) 



with r) = — 1. In the above formula, <Ji are Pauli op- 
erators. This representation allows us to calculate in a 



straightforward way, using trace properties of Pauli oper- 
ators, the joint probability of measuring the two photons 
in respective polarizations u a and Uf,: 

p(u a , Uf,) = jTr[|3(i + u a • &) ® (i + u 6 ■ <x)] 

= -[l + 77u„ -u 6 ] (3) 

where <r denotes a three-component vector composed of 
the Pauli operators. Thus the joint probability for po- 
larization settings u a and Uf, is a function of the an- 
gle a between the corresponding Bloch vectors, given by 
cos a = u a • Uft. When performing a common transfor- 
mation on both measured polarizations u Q and Uf,, the 
probability is left unchanged. It is seen that the visibil- 
ity of polarization correlations is given by \tj\, reaching 
100% for the singlet state. The density operator defined 
in Eq. © is positive definite for — 1 < r\ < |. Thus 
although perfect anticorrelations in the polarization de- 
gree of freedom are possible, there is no physical quan- 
tum state in which the two photons would always have 
identical polarizations in an arbitrary basis. 



III. SEMICLASSICAL BOUNDS 

Let us consider a pair of light beams distributed be- 
tween two parties a and b. Because of the employed de- 
tection scheme, it will be sufficient to take into account 
only the intensity and the polarization degrees of free- 
dom. Classically, a general state of polarization is de- 
scribed by a four-component Stokes vector. For a fully 
polarized beam the component corresponding to the to- 
tal intensity is unambiguously defined by the remaining 
three coordinates. In order to retain correspondence with 
the quantum mechanical case, we will order these coordi- 
nates identically as in the Bloch vectors and we will call 
here the resulting three-component objects generalized 
Bloch vectors. 

Thus an individual preparation of the two beams is 
described by generalized Bloch vectors s Q and S&, whose 
length corresponds to the intensity, and the direction — to 
the polarization state of the respective beam. If the com- 
plete knowledge about the beams is unavailable, the joint 
state is described by a probability distribution P(s a ,Sf,). 
In the following, we shall assume that the state is in- 
variant with respect to an arbitrary common polariza- 
tion transformation P(fts a ,ftsb) = P(s ,S(,), where ft 
belongs to the SO (3) group of proper rotations in three 
real dimensions. 

The measurement is performed by selecting polariza- 
tion components corresponding to u Q and u/, in the re- 
spective beams, and detecting clicks using single pho- 
ton counting modules. Assuming that the detector fires 
whenever one or more photons are detected, the click 
probability for the classical light of normalized inte nsity 
s entering the detector reads n(s) = 1 — exp(— s) [22j |. 
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Therefore, the probability of a coincidence event in the 
semiclassical model is given by 



p(u a ,u b ) = 

ds a J ds fc P(s a ,s f ,)n(i(|s Q |+u a -s Q ))n(i(|s 6 |+u 6 -s h )) 



(4) 



Detector losses are described by multiplying the argu- 
ment of IT(s) by the detection efficiency, which can be 
included in the above formula by appropriate rescal- 
ing of the probability P(s a ,Sb). Let us now represent 
s a = s a l~2e, where e is a reference unit vector, 11 is a 
certain rotation, and s a — \s a \. Writing the integration 
measure as ds a = ds a dtt and making use of the assumed 
rotational invariance of P(s a ,s b ) allows us to write 

p(u Q ,u b ) = J ds a J dfl J ds b P(s a e, n _1 s b ) 

x n(| So (i + u a • n e ))n(i(| S6 | + u b • S6 )) . (5) 

Substituting s b — * l~2s b yields: 

p(u a ,u b ) = J ds a J ds b P(s a e 7 s b ) 

J dOn(is Q (l + u Q -^e))n(i(|s b |+u 6 -ns b )). (6) 

It will be convenient to assume in the following that the 
measure over 17 is normalized to one, J dfl = 1. Us- 
ing the explicit form of the count probability n(s), the 
integral over 17 in Eq. ^ can be simplified to 

J dl7n(±s Q (l + u Q • fie))ll(i(|s 6 | + u b • 17s b )) = 
!_ _L(i _ e -KI) _ JL(i _ e -M) + e -(l-.l+l- t |)/3 Z / j 



S6 



(7) 



where the only term depending on the polarizations u a 
and u b being measured reads 



U= I dfle~'- SaUa ' ne+Ub ' nSb]/2 . 



(8) 



This integral can be calculated explicitly when s a = — s b , 
leading to a closed expression for the coincidence proba- 
bility p(u a , u b ) of the form 



P(u a ,u b ) 



1 



sinh(ssin§) , 
2sm ^ 



, (9) 



where a is the angle between the vectors u and u b and 
s = \s a \ = |s b |. In Fig. Q] we depict the dependence 
of the p(u a , u b ) on the angle a for several values of the 
intensity s. In order to aid comparison, the graphs have 
been rescaled by the the probability p tot of a coincidence 
event when the polarizers are removed from the setup. It 
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FIG. 1: The coincidence probability as a function of the an- 
gle a between the Bloch vectors of measured polarizations 
for the singlet state (dotted, black) and a rotationally invari- 
ant classical state of two beams with perfectly anticorrelated 
polarizations and identical intensities, equal to s = 0.0025 
(dashed-dotted, light gray), s = 1 (dashed, dark gray), and 
s — 10 (solid, black). The curves are normalized by the prob- 
ability ptot of a coincidence event when entire beams without 
selecting polarization components are detected. 




FIG. 2: The visibility of polarization correlations in the 
semiclassical model when both beams have equal intensities 
\s a \ = \sb\ — s and either identical or orthogonal polariza- 
tions. 



is seen than the correlations are much weaker than in the 
quantum mechanical case of a singlet state, and for larger 
intensities they loose the purely harmonic form. Note 
that choosing s a = s b in Eq. ([5]) leads to an expression 
analogous to Eq. ^ with sin ^ replaced by cos |f. 

We will now derive bounds on IA for arbitrary general- 
ized Bloch vectors s a and s b , which will be used to an- 
alyze the visibility of polarization correlations in a more 
general case. In order to find a lower bound on U let 
us start from an observation that the integral in Eq. ([8]) 
is invariant under the substitution 17 — ► Olio- Taking 
S7o to be a it rotation around the axis perpendicular to 
the plane spanned by e and u b implies that the value 



FIG. 3: A bound on the visibility of polarization correlations 
achievable using classical light, depicted as a function of the 
beam intensities |s a | and |sb|. The bound is tight for |s a | = 
1 st, 1 - In the white regions near axes the bound exceeds the 
value |, conjectured to be the general visibility limit in the 
semiclassical regime. 



FIG. 4: The path drawn by the Bloch vector of a linearly 
polarized beam sent through a pair of counter-rotating wave- 
plates with phase delays 7r and arccos ^= . This transforma- 
tion applied to two beams with initially orthogonal polariza- 
tions produces an ensemble equivalent in the limit of very 
weak intensities to a rotationally invariant one. 



of hi does not change under simultaneous substitutions 
e — ► — e and u/, — > — uj,. This allows us to write: 

If _ I J rf ^ e -(s a u a -f2e+u b -f2s b )/2 + e (s a u a -ne+u b -r2s b )/2j 

= J dfl cosh(|s a u Q ■ fie + u 6 ■ fls b \/2). (10) 

Using the convexity of the hyperbolic cosine function and 
the inverse triangle inequality we obtain 



U > cosh / dtt\s a u a ■ He + u b ■ S~2s;,|/2 



> cosh (^J dU{\s a u a ■ fie| - |u b • fts b |)/2^ . (11) 

Evaluating the two integrals over $7 yields 

U > cosh[(|s Q | - |s 6 |)/4]. (12) 

An upper bound on U can be found with the help of the 
Cauchy-Schwarz inequality 



U < 





sinh |s a | sinh |s&| 
|s a ||sb| 



(13) 



The two derived bounds allow us to estimate from above 
and from below the coincidence probability p(u a , u b ) over 
all possible settings of u a and u;,. 



IV. CORRELATIONS VISIBILITY 

Let us first discuss the limit of weak intensities 
|sfc|,|sb| <C 1, when the probability of a click on a de- 
tector can be approximated by n(s) w s. In this case, 
the probability of a coincidence event can be written as 



p{u a ,u b ) 



1 



(|s a ||s 6 |> + -(s a -s fe >cosa (14) 



where the angular brackets denote the average (...) = 
J ds a J dsi,P(s a , s;,) . . .. The Cauchy-Schwarz inequality 
means that |(s a -S(,)| < ^|s a ||sf,|), which implies that the 
visibility cannot exceed |. 

In a symmetric scenario, when both the beams always 
have orthogonal polarizations and equal intensities |s a | = 
|s(,| = s , we can use Eq. © to calculate the visibility of 
polarization correlations 



V 



sinh s 



s + (2s - 3) sinh s + 2(s - 2) cosh s + 4 ' 



(15) 



This quantity, shown in Fig. [21 is a monotonically de- 
creasing function of the intensity s, with the maximum 
value V = | reached in the limit s — > 0. Note that in 
this scenario the bounds given in Eqs. (fi"2|) and (fTS"]) are 
tight. 

More generally, we can use the inequalities derived in 
Eqs. (|T2]) and (fl3)) to find an upper bound on the visibility 
of polarization correlations for arbitrary intensities |s Q | 
and |s;,|. This bound is depicted in Fig. [31 It is seen that 
the visibility indeed does not exceed i except regions 
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FIG. 5: The experimental setup. NDF, a set of neutral den- 
sity niters; A/2, A/4, half and quarter wave plates; WP, a wave 
plate introducing an arccos rad phase shift; PBS, polariz- 
ing beam splitters; SPCM, single photon counting modules. 
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FIG. 7: The coincidence count rates as a function of the ro- 
tation angle <p of the polarizer PBS2 placed in the path of 
the beam b. The three panels display results for scans of the 
polarization of the beam 6 along the great circles in the xz 
plane (top), the yz plane (middle), and the xy plane (bot- 
tom). The solid lines depict least-squares sinusoidal fits. The 
insets specify the polarization of the beam a for each of the 
datasets along with visibilities obtained from the fits. 



FIG. 6: The three great circles on the Bloch sphere specifying 
the polarization scanned for the beam b to verify polarization 
correlations. The blue circle in the xz plane was scanned for 
the horizontal <-> and the vertical | polarizations of the beam 
a. The red circle in the yz plane was scanned for circular 
polarizations O and O of the beam a. The green circle in the 
xy plane was scanned for diagonal polarizations f and \ of 
the beam a. 



when one of the intensities |s a | or |sb| is small. Let us 
have a closer look at the limit when one of the intensities 
approaches zero. For concreteness, let us take |s a | — > 0. 
It is then justified to expand 

J <m(l-\u a -ns a )e-" b ns »' 2 . (16) 

This expression can be evaluated analytically to yield 



2 IsJ 1 , 

U « ] — | smli — + - — ^(Ua • u b )(s a e • s b ) 

Sfc 2 \ s b\ 



x [ cosh — — r sinh 

V 2 |s 6 | 2 



(17) 



product (u a • U(,)(s a e • s&) = ±|s a ||s&|. It can be veri- 
fied by a straightforward yet lengthy calculation that the 
visibility of polarization correlations obtained using this 
expression approaches | when |s a | — * 0. This shows that 
the upper bound derived with the help of the inequalities 
given in Eqs. (|12|) and (I13|) is not tight and suggests that 
the limit ^ might be universal for arbitrary intensities 
|s a | and |s b |. 

Let us note that the bound ^ does not discriminate 
between positive and negative polarization correlations. 
The quantum mechanical counterpart a statistical ensem- 
ble discussed above is a situation when each of the two 
photons has a well defined polarization state, and only 
statistical averaging is permitted, i.e. the pair is prepared 
in a separable state. For Werner states introduced in 
Eq. ([2]) the separability criterion takes a form symmetric 
with respect to the sign of the parameter r\: —\<ri<\ > 
[HI, also leading to the visibility limit of -| for polariza- 
tion correlations. 



EXPERIMENT 



It is seen that the minimum and the maximum val- 
ues of the coincidence probability are obtained when the 



The semiclassical limit of polarization correlations can 
be demonstrated experimentally using a pair of orthogo- 
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nally polarized light beams, submitted to identical polar- 
ization transformations selected randomly according to 
the Haar measure on the rotation group SO (3). Mechani- 
cal implementation .2_4j of such depolarization would lead 
to excessively long averaging times. However, for low in- 
tensities, when the probability of a count is proportional 
to the incident intensity, the polarization measurement 
is sensitive only to certain moments of the distribution 
of output polarizations: 

p(u a ,u b ) | ((|s a ||s a |) +u a • (|s a |sfc) +u b • (s a |s b |) 

+u a • (s^s h )u 6 ) (18) 

These moments can be reproduced by a discrete set of 
Bloch vectors, or a continuous one-parameter family. The 
latter solution is a convenient choice for an implementa- 
tion in an optical setup. An exemplary realization is a set 
of two waveplates with phase delays equal respectively to 
arccos rad and tt rad that rotate in opposite directions 
with the same constant angular velocities. For a suitably 
selected input linear polarization the output Bloch vector 
draws a curve depicted in Fig. |4j which yields the output 
ensemble with desired statistical properties. 

The complete experimental setup is presented in Fig. El 
The light source was a highly attenuated pulsed beam 
from a mode-locked Tksapphire oscillator with the cen- 
tral wavelength 775 nm and the full width at half- 
maximum bandwidth 10 nm, corresponding to subpi- 
cosecond pulse duration. The polarizing beam splitter 
PBS1 followed by a silver mirror was used to prepare 
two parallel beams in horizontal <-*■ and vertical | polar- 
izations with respect to the plane of the setup. The half 
wave plate inserted before PBS1 served to equalize the 
intensities of the two beams. The two wave plates intro- 
ducing depolarization were mounted on rotation stages 
driven by stepper motors. In order to produce the correct 
output ensemble, the axes of the waveplates were initially 
aligned parallel to each other. For each measurement 
point, the waveplates were rotated in opposite directions 
with constant angular velocities, completing an integer 
number of 10 full rotations. The two emerging beams 
were transmitted through separate quarter wave plates 
A/4 followed by polarizers, which selected individual po- 
larization components, and then focused on the active 
areas of single photon counting modules SPCM (model 
Perkin Elmer SPCM-AQ-131). The input beam was at- 
tenuated to an intensity at which single count rates on 
the SPCMs with removed polarizers were approximately 
equal to 2 x 10 5 Hz. Given the laser repetition rate of 
80 MHz this yields the intensity parameter s w 2.5 x 10~ 3 . 

In order to verify the correlations, we chose six differ- 
ent polarizations for the beam a: horizontal <->, vertical 
I, diagonal at ±45 degrees y and \, as well as right- and 
left-circular O and 0. For each of these six polarizations, 
we measured coincidence count rates as a function of the 
polarization selected for the beam b sweeping one of three 
great circles on the Bloch sphere, shown in Fig. [5] The 
great circles were scanned by setting the quarter wave 



plate in the beam b to an the appropriate orientation 
or removing it altogether for the scan of the xz circle, 
and rotating the polarizer PBS2 in 10° steps, which cor- 
responds to 20° increments on a great circle. For each 
setting, photocounts were collected over a 60 s time in- 
terval, and both single and coincidence count numbers 
were recorded. 

The measured polarization correlations are depicted in 
Fig. The data sets were fitted using the least-squares 
method with sinusoidal patterns to obtain corresponding 
visibilities, specified in the graph insets. The experimen- 
tally determined values are close to the predicted value 
of -| . Small deviations can be attributed to non- uniform 
depolarization of the two beams. This imperfection could 
also be noticed as weak dependence of single count rates 
on the measured polarization, with relative variations up 
to 1.5%. 



VI. CONCLUSIONS 

We investigated the strength of polarization correla- 
tions that can be achieved using classical optical fields 
and detectors without photon number resolution. We 
showed under certain symmetry assumptions that there 
is a substantial gap between the semiclassical regime and 
the correlations allowed by quantum mechanics. It is 
worthwhile to note that the maximum value g of the vis- 
ibility of polarization correlations is lower than the bound 
of -4= in local realistic theories required to satisfy stan- 
dard Bell's inequalities. This illustrates that the semi- 
classical description of optical fields is only a particular 
case in the general class of local realistic theories. 

From a more general perspective, the gap between 
classical and quantum theories calls for a deeper under- 
standing of the role of quantumness in information pro- 
tocols based on quantum interference and entanglement. 
This research direction has been shaped by Krzysztof 
Wodkiewicz in his last works [25|, [26|, H3, One in- 
triguing problem is to what degree and under what as- 
sumptions classical theories can describe quantum proto- 
cols for information processing, especially when one uses 
input states that have classical counterparts. This issue 
has been analyzed thoroughly for the case of continuous- 
variable quantum teleportation [25| . which motivates 
looking into other scenarios, one interesting candidate 
being quantum cloning of spin-coherent states 26]. Fur- 
thermore, despite dramatic differences between predic- 
tions of quantum mechanics and classical theories, it is 
possible to draw insightful parallels between classical and 
quantum interference phenomena (27| . leading to feasi- 
ble experimental observations in the optical domain [28| . 
This opens up a question whether one can find ways to 
employ classical interference in some applications that 
are currently thought to depend indispensably on quan- 
tum resources. 
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